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ABSTRACT 

We investigate the interaction of nonthermal ions (protons and nuclei) accelerated in an ultrarelativistic blast- 
wave with the pre-existing magnetic field of the medium into which the blastwave propagates. While particle 
acceleration processes such as diffusive shock acceleration can accelerate ions and electrons, the accelerated 
electrons suffer larger radiative losses. Under certain conditions, the ions can attain higher energies and reach 
farther ahead of the shock than the electrons, and so the nonthermal particles can be partially charge-separated. 
To compensate for the charge separation, the upstream plasma develops a return current, which, as it flows 
across the magnetic field, drives transverse acceleration of the upstream plasma and a growth of density contrast 
in the shock upstream. If the density contrast is strong by the time the fluid is shocked, vorticity is generated at 
the shock transition. The resulting turbulence can amplify the post-shock magnetic field to the levels inferred 
from gamma-ray burst afterglow spectra and light curves. Therefore, since the upstream inhomogeneities are 
induced by the ions accelerated in the shock, they are generic even if the blastwave propagates into a medium of 
uniform density. We speculate about the global structure of the shock precursor, and delineate several distinct 
physical regimes that are classified by an increasing distance from the shock and, correspondingly, a decreasing 
density of nonthermal particles that reach that distance. 

Subject headings: acceleration of particles — cosmic rays — gamma-rays: bursts — plasmas — shock waves 



1. INTRODUCTION 

The standard model for gamma-ray burst (GRB) after- 
glow emission requires nonthermal acceleration of electrons 
to ultrarelativistic energies. This can take place in the "ex- 
ternal" relativistic shock wave driv en into the c i rcumburst 
mediu m by the burst ejecta (e.g., iPiranl l2005at iMeszaroi 
2006, and references therein). Candidate mechanisms for 
shock acceleration, such as the first order Fermi mechanism — 
also known as diffusive shock acceleration (PSA; e.g.. 



Beii[T97a iBlandford & Ostrikeri[l978T: felandford & Eichlen 
1987) — can accelerate all charged particles, including the 
positively charged ions, for whose acceleration in GRB 
afterglows we still do not have direct observational evi- 
dence. While theoretical models for particle acceleration 
in shocks have largely been based on analytic or semian- 
alytic prescriptions for particle injection into the accelera- 
tion cycle and its scattering in the random magnetic field 
of the shock (e.g., [A chterber g etal]|2001t lEllison & Double! 
2002; Niemiec & Ostrowski 2006, and references therein), re- 
cent two-dimensional first-principles particle-in-cell simula- 
tions of unmagnetized pair shocks have delivered the first en- 
couraging direct evidence fo r the formation of the signature 
nonthe rmal power-law tail dKeshet et al.l l2008t ISpitkovskvl 
2008b). The power law tail seems to arise in tandem with a 
self-sustained electromagnetic structure in the shock precur- 
sor, where the nonthermal particles that move ahead of the 
shock excite plasma instabilities. In their nonlinear devel- 
opment, the instabilities generate a magnetic field that scat- 
ters the nonthermal particles back toward the shock and thus 
closes the acceleration cycle. Our aim here is to explore the 
consequences of the formation of a power-law tail of accel- 
erated particles for shock hydrodynamics when a weak, large 
scale magnetic field is present in the shock upstream. 

1 The University of Texas, 1 University Station C1400, Austin, TX 78712 

2 Theoretical Astrophysics, Mail Code 130-33, California Institute of 
Technology, 1200 East California Boulevard, Pasadena, CA 91 125. 



Inverse Compton cooling of the nonthermal electrons im- 
poses an upper limit on the energy of electro ns that are ac- 
celerated by DSA (e.g.. iLi & Waxmanll2006l) . Since higher 
energy ions reach farther from the shock into the shock up- 
stream, there may exist a region in the shock upstream devoid 
of nonthermal electrons that still contains nonthermal ions; 
we restrict our analysis to this region. In the reference frame 
in which the shock upstream is at rest, the nonthermal ions are 
beamed forward within the narrow angle ~ T" 1 , where r> 1 
is the shock Lorentz factor. The upstream plasma responds 
to the forward-beame d nonthermal ions by developing a re- 
turn c urrent (see, e.g.. lThode & Sudanlll975l : ISpicer & Sudani 
1984, and references therein) which flows toward the shock. 
The return current cancels the charge separation between the 
nonthermal ions that reach farther from the shock and the non- 
thermal electrons that are confined closer to the shock by the 
inverse Compton losses. If the upstream plasma is weakly 
magnetized, the flow of return current across a weak mag- 
netic field accelerates the upstream plasma transversally, i.e., 
perpendicular to the direction of shock propagation, which 
creates inhomogeneities in the plasma dReville et al.l 120061 
iPelletier. Lemoine. & Marcowithll2008l); this is a relat ivistic 
version of the mechanism discussed by Bell (20041 12005b . 

We show that even when the upstream magnetic field 
is weak, these density inhomogeneities may be sufficient 
to generate a stronger magnetic field in the shock down- 
stream via vorticity production at the shock transitio n 
dGoodman & MacFadvenl 120071: ISironi & Goodmanl |2007|) . 
This process converts the bulk kinetic energy of the shock 
downstream into vortical and magnetic energy. An advan- 
tage of this process is that it can take place even when 
the circumburst medium is uniform and therefore it is 
relevant for, as required by observations, both long and 
short GRBs (Nakar 2007). Alternative sources of vortic- 
ity that were previously discussed in the literature are pre- 
existing density inhomogeneities in situations where an ul- 
trarelativistic blastwave propagates into a clumpy stellar 
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wind ( ISironi & Good man 200 7|), or an initially anisotropic 
blastwave (Milosavljevic, Nakar, & Zhang 2007). Nonther- 
mal particle streaming may generate or amplify a mag- 
netic field in the shock u pstream dBell & Lucekl 120011: 
Milosavljevic & Nakar 2006; A. Spitkovsky, pri v. comm., but 
also see lPelletier. Lemoine. & Marcowithll2008h . The relative 
efficiencies of the direct, nonthermal particle-driven field am- 
plification in the shock upstream, and the amplification in the 
shock transition via vorticity generation in the presence of 
streaming-generated density inhomogeneity will depend on 
how the shock energy is partitioned between the thermal fluid 
and the nonthermal accelerated particles and on the efficiency 
with which magnetic field is generated by cosmic ray-driven 
turbulence in the shock upstream. 

Throughout this work we assume that in ultraretivistic col- 
lisionless shocks, DSA accelerates ions to high nonthermal 
energies. While there is direct observational evidence of elec- 
tron acceleration in such shocks, the assumption that pro- 
tons are accelerated also may not necessarily be justified and 
deserves further investigation. Specifically, on theoretical 
grounds, particle acceleration should be inefficient when a 
compressed oblique large scale upstream magnetic field dom- 
inates particle dynamics in the shock downstream. In this 
case shock compression leads to a downstream magnetic field 
that is almost perpendicular to the shock normal. Then, a 
particle with a gyroradius smaller than the coherence length 
of the field cannot complete two full DSA cycles across the 
shock before becoming irrever sibly bound to and e scaping 
into the shock downstream (e.g., Lemoi ne et al.l2006h . There- 
fore, a field that contains large-amplitude magnetic field per- 
turbations on small scales seems to be required for accelera- 
tion (e.g.,|O strowski & Bednarz 2002; Niem iec & Ostrowskil 
12004 l2006t iNiemiec et al.1 120061) . Even in the conditions in 
which a shock propagates in a weakly magnetized medium 
and generates magnetic turbulence that serves as a scattering 
agent for energetic particles, it cannot be taken for granted 
that the resulting particle spectrum will extend up to the high 
Lorentz factors that are required for the mechanism that we 
present. These potential limitations should be born in mind in 
what follows. 

This work is organized as follows. In § 12.11 we analyze 
charge separation in the shock precursor and relate the non- 
thermal ion density ahead of the shock to the maximum en- 
ergy to which nonthermal ions are accelerated in the shock. 
In § 12.21 we study the motion of the upstream medium in re- 
sponse to the Lorentz force associated with the nonthermal 
ion streaming across a weak magnetic field. In § 12. 31 we de- 
termine the conditions under which nonlinear density inho- 
mogeneities can result and estimate the energy density of the 
magnetic field produced by vortical motions that are excited 
a s the shock transition sweeps past the inhomogeneities. In 
§ 12.41 we examine the production of a small scale reversing 
magnetic field structure that was assumed in § 12.21 in the 
weak large-scale magnetic field of the circumburst me dium 
by the nonr esonant streaming instability introduced by iBelll 
(2004, 2005). In §[3] we speculate about the global structure 
of the shock precursor, and delineate several distinct physical 
regimes that are classified by an increasing distance from the 
shock. In §2] we summarize our main conclusions. 

2. INTERACTION OF THE SHOCK PRECURSOR WITH 
THE UPSTREAM PLASMA 

2.1. Charge Separation 



We consider a relativistic blast wave described by Lorentz 
factor r in which electrons and positive ions are accelerated 
by DSA. Unless otherwise explicitly noted, all quantities are 
defined in the rest frame of the upstream. The maximum 
energy to which the electrons are accelerated, 7 e . ma x'«pC 2 , is 
limited by inverse Compton cooling as the ele ctrons traverse 
the sh ock upstream during a DSA cycle (e.g.. iLi & Waxmanl 
2006). The ions experience much smaller radiative losses 
and can be accelerated to higher energies, ji, max niiC 2 , than 
the electrons. Here, m, is the ion mass, and j e and 7, de- 
note the electron and ion Lorentz factors, respectively. Hence, 
the accelerated ions can travel farther ahead of the shock 
than the electrons. If 7,, ma x'W/ > 7e.max'«?, the shock pre- 
cursor consisting of accelerated particles will not be charge 
neutral far from the shock, at distances that the nonthermal 
electrons cannot reach. There, charge neutrality is recovered 
by a small polarization of the upstream, unshocked plasma. 
The polarization is established by th e flow of a return cur- 
rent in the upstream p lasma (see, e.g.. lThode & Sudanll 1975b 
Spicer & Sudan iT984L and references therein). This return 
current drives the growth of upstream plasma inhomogeneity 
dBellll2 004, 2005} that is the subject of this work. 

The maximum distance A(7) that a nontermal particle with 
Lorentz factor 7 can reach from the shock is determined by 
the magnetic field structure in the shock upstream. The mag- 
netic field can come from the weak pre-existing magneti- 
zation of the circumburst medium, or can be generated by 
plasma self-organization resulting from an interaction with 
the streaming of nonthermal particles. Depending on the 
structure of the magnetic field, nonthermal particle orbits are 
either deflected in a circular fashion if the magnetic field is 
coherent, or execute random walks if the field is tangled, 
i.e., reverses on scales A smaller than the gyroradius of the 
nonthermal particles divided by the shock Lorentz factor, 
r g /T ~ 7«c • /eBT, where B is the RMS field strength. Since 
the magnetic field can have different strength and structure on 
different scales B(X), we employ the subscript "def" to de- 
note the magnetic field length scale that produces the largest 
(circular or random) deflection of the nonthermal particle. 

When a nonthermal particle is deflected by angle > T" 1 
from the shock normal, its velocity parallel to the direc- 
tion of shock propagation becomes smaller than that of the 
shoc k and the shock can catch up with the particle (see, 
e.g.. lAchterberg et al.ll200ll and references therein). A par- 
ticle that has been overtaken by the shock can be scattered 
back into the shock upstream, as is necessary to close the 
DSA cycle. Therefore, working in the reference frame in 
which the shock upstream is at rest, for the coherent field we 
have A(7) ~ r g /T 3 , whereas for the tangled field we obtain 
A(7) ~ fg/Adefr 4 , where r g = -fmc 2 / 'eBdef- Since we do not 
estimate the magnetic field strength and geometry from first 
principles, we parametrize the effect of the deflecting, possi- 
bly small scale magnetic field by the maximum Lorentz factor 
of a nonthermal ion 7 p , ma x, which we here take to be a proton, 
that can be accelerated with the aid of the magnetic field. We 
separately treat the cases in which the deflecting field is coher- 
ent and tangled. In the case of a tangled field, we assume that 
the field is tangled on spatial scales Ad e f <C r„/T that are in- 
dependent of the distance from the sh ock A. This assumption 
may not be valid in ge neral (see, e.g., Milosavljevic & Nakarl 
l2006tlKatzetai1l2007h . in which case our analysis should be 
interpreted as applying to the highest energy accelerated parti- 
cles, and to the magnetic field component, with its own char- 
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acteristic correlation length Adef, that has the greatest influ- 
ence on the dynamics of these particles. 

We also assume that the blastwave propagates into a quasi- 
homogeneous medium with average density p ~ m p n. The 
highest energy protons can reach the farthest from the shock 
to a distance A p (-f p _ m . dx ) ~ R/&T 2 , where R is the radius of 
the blastwave, and the factor of one eighth is peculiar to adi- 
abatic spherical ultrarela tivistic blastwaves propa gating into 
uniform density media (Blandford & McKee' 119761) . Let U s h 
denote the energy density in radiation at radius R in the shock 
frame which is a fraction e ra d = (47r/3)/? 3 t/ s h/-£'tot of the to- 
tal isotropic equivalent energy of the blastwave E tot . Equat- 
ing the inverse Compton power |0rc 7 2 sh t/ s h, where 7f ,. s h ~ 
j e /T is the electron Lorentz factor in the shock frame, ctt 
is the Thompson cro ss section, and we n eglected the Klein- 
Nishina effects (e.g. jLi & Wax man 2006), to the energy gain 
~ 7e,sh»JeC 2 per DSA cycle of duration ~ A e:Sh /c ~ TA e /c, 
we obtain the maximum Lorentz factor to which electrons 
can be accelerated, in spite of the cooling, as a function of 
the Lorentz factor to which the non-cooling protons can be 
accelerated 



7e,n 



4.4c'/3 m y 2 7 ;/i x r'/3 



1/2. 
£ rad f 



V3P 



1/6 1/2 



(coherent field) , 



2.i cy^JVJLw 



£ l/3 
' rad 



P j££ 

1/3 1/3 



(tangled field) , (1) 



which, as usual, is expressed in the rest frame of the shock 
upstream. Here, we have eliminated the blastwave radius R 
via the relation R = (17£ , tot /87iT 2 /5c 2 ) 1/ ' 3 applicable to an adi- 
abatic spherical blastwaves propag ating into uniform density 
media (Blandford & Mc Ked 19761) The magnetic field length 
scale Adef (in the case of a tangled field) and strength squared 
B\ e{ do not appear in equation ((TJ because they have been ex- 
pressed in terms of the maximum proton Lorentz factor that 
can be accelerated by deflection in a field with these proper- 
ties via Bdef = lhp,maxtn p c 2 /e)/r 3 ]/(R/8r 2 ) (coherent field) 
and AdefflL = [(>,max'« / ,c 2 /e) 2 /r 4 ]/(i?/8r 2 ) (tangled field). 

The nondimensional distance from the shock that the most 
energetic electrons can reach, relative to the distance that the 
protons can reach, can be expressed as 

A c (7e,max) 



■^cool ; 



Ap(7/?,max) 

4.4 c'/VT 1 / 3 
-v 1/2 e 1/2 o 1 / 3 £ 1/6 m I/2 ( j 1/2 

7 c 4 /V /3 r 4 / 9 

2/3 2/3_ 4/9c ,2/9 2/3 2/3 



(coherent field), 



(tangled field) . (2) 



For typical values of e rad , Etot) T, and «, this becomes 



-^cool ' 



V 460 / \0.l) 



l / 2 /erad\-V2 / £ tc 



100 



1/3 



1 cm 



V 10 53 erg 

-1/3 



-1/6 



/ 7p,max \ 

V 440 / 



(coherent field), 

-2/9 



7p,max V 2/3 (^±Y 2/1 ( Exot 

V0.1/ V 1053 er S, 

4 / 9 , - , -4/9 

(tangled field), (3) 



-Y 

100/ 



1 cm 



which shows that for 7 p . max 3> 10 3 , regardless of the magnetic 
field geometry, the electrons are unable to travel as far from 
the shock as the protons and x coo i <C 1 in either case. Note that 
for a coherent confining field 



7p,max~7 X 10 5 



r 



B, 



def 



( E u 



10"V V 10 53 erg 



1/3 



100 



1/3 



1 cm 



1/3 



(coherent field), (4) 



implying that if the confinement is dominated by a pre- 
existing microgauss upstream magnetic field we expect that 
7^.max ~ 10 5 - 10 6 , while a larger value is expected if the 
shock precursor generates a magnetic field in which the tan- 
gled component dominates particle deflection. Therefore we 
always expect that jt coo i < 1. In what follows, we ignore 
the interior region A < A e (j e , max ) populated by nonthermal 
particles of either sign and assume that the exterior region 
A > Aeije , max ) contains only nonthermal ions, which we have 
assumed to be protons. 

The nonthermal particle concentration at a given distance 
ahead of the shock will be dominated by the lowest energy 
particles that reach that distance. Thus, we can relate the non- 
thermal proton density n ntp in the shock upstream to their en- 
ergy distribution 



Mltp 



(7)' 



1 



4tt/? 2 A(7) dlnj 1 



(5) 



where JVjj? (7) is the total number of nonthermal protons in 
the shock upstream with Lorentz factors less than 7. 

We assume that the accelerated particle spectrum in the 
shock downstream resulting from DSA contains equal energy 
on all energy scales, dN^ /d-r ex 7 2 dLemoine & Pelletierl 
120031; lEllison & Doublell2002l 12004 cf. e.g.,); this assump- 
tion is not critical as somewhat steeper spectra, expected 
when particle scattering in the shock downstream is not 
isotro pic (e.g.. iKeshet & Waxman| [2005; Lemoine & Revenul 
120061) . lead to similar conclusions. Ignoring the energy in 
nonthermal protons located in the upstream, we may normal- 
ize the spectrum by requiring that 



7p,» 



T/'.min 



2«iv ntp 

jm p c dj- 

«7 



(6) 



where e nt is the fraction of the total energy in the accelerated 
protons. The minimum Lorentz factor of the nonthermal par- 
ticles, resul ting from a single scattering across the shock is 
7p,min ~ T 2 dGallant & AchterbergTi 999). but the dependence 
on j p>m in is only logarithmic. 

To relate A^p P) to A^f p own) , let P ~ \ denote the probability 
that a particle in the shock upstream, upon returning to the 
shock, is scattered again into the shock upstream. The down- 
stream then contains the particles that are not re-scattered into 
the upstream, but are entrapped within and are advecting with 
the downstream. Conservation of nonthermal proton flux at 
the shock transition then implies 3 



dN 



(down) 



ntp 



d-p) 



tdyn dN, 



(up) 
ntp 



face (7) rfln 7 



3 This corrects the inadequate assumption A^"p'(7) 
had made in Milosavljevic & Nakar (2006). 



(7) 



(7) that we 
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Here f acc is the time between the emission of a nonthermal 
proton from the shock and its reabsorption in the shock 



4cc (7) ' 



(i-%h) 1/2 r 2 A( 7 ) 

/%8h C 



(8) 



where /3u sh ~ 5 is the average component parallel to the di- 
rection of shock propagation of the velocity of nonthermal 
protons in the shock frame, and fdyn ~ R/c is the dynamical 
age of the blastwave. From this, 



dN 



(up) 



ntp 



1 



(down) 
ntp 



(9) 



din 7 \7p,maxy d\n-y 

where s = a! In A/a! In 7, so that 5=1 for the coherent field, 
and s = 2 for the tangled field. This shows that the integral 
spectrum of all the nonthermal particles located in the shock 
upstream is substantially harder than that in the downstream. 

Substituting equation © in equation (0, we find that the 
total density of the accelerated protons equals 



'ntp 



o.04 ent p 2 / 3 4{ 3 r 4 / 3 



Ara„c 2 / 3 



7w 



.A 



(coherent field) , 



0.1e nt p 5 /^f r 8/3 



.1/3,. 



(tangled field) , (10) 



where again A(7) = (7/7p.max) s ^/8r 2 and we have defined 
A = ln( 7p , max / > . min ) - 5-20. Let x = A/(/?/8r 2 ) denote 
the nondimensionalized distance from the shock; in terms of 
this parameter, 

, , 0.3 /-,,,, ^ 4 



r 

Too 



) 4 (t=*) (&)•<"> 



where q = s~ , so that q = 1 for the coherent field and q =_= 
for the tangled field, and we have set A = 10. Equation dl jj 
gives the nonthermal proton density in the rest frame of the 
shock upstream. For large V or small 7 / >, max , this density can 
formally exceed the upstream density; in reality, n ntp may be 
limited by kinetic plasma instabilities. In the rest frame of the 
nonthermal proton shock precursor, the density ratio n ntp /fi is 
reduced by factor T" 2 . 

2.2. Dynamics 

Here we study the response of the shock upstream to the 
streaming of the nonthermal protons in the shock precursor. 
The streaming of nonthermal protons with density « ntp will 
give rise to a return current 



Jret(A) ■ 



'ntp 



(A)R, 



(12) 



where R is a unit vector in the direction of propagation of 
the shock. The return current cancels the charge separation 
among the nonthermal particles resulting from the inverse 
Compton cooling suppression of the acceleration of nonther- 
mal electrons. Consider a fluid element in the upstream 
medium containing a magnetic field B that reverses direction 
on length scales A. This may not be the same magnetic field 
scale that is the most effe ctive at deflecting the orbits of the 
nonthermal particles, see § 12.11 A magnetic field that reverses 
on the relevant scales may exist in the upstream at the outset. 
If not, as we show in § 12.41 it can be produced by the return 
current from an initially coherent microgauss field. Let [i de- 
note the sine of the angle between the magnetic field and the 
direction of shock propagation. 



As the return current flows across the field, the plasma expe- 
riences a Lorentz force acting perpendicular to both the direc- 
tion of shock propagation and the field direction. Let r denote 
the time measured since the arrival of the photon shell; in this 
time coordinate, the shock transition sweeps a fluid element 
at t = R/8T 2 c. The equation of motion for the transverse La- 
grangian displacement of the fluid element in response to the 
Lorentz force, valid while the motion is Newtonian, i.e., while 
\d£/dr \ <c, reads 



dr 2 



, Jret(r) x B 

pc 



(13) 



From this, assuming that B ~ B sin(27rr • B / A)B, where r is 
a spatial coordinate and B is peak field, the amplitude of the 
transverse acceleration is 



d\ 
dr 1 



1.7 ttecejfp 1 !^ 



vA 



(14) 



where tj> = 2irr ■ B/A, and €b = B 2 1%-npc 2 is the energy den- 
sity in the magnetic field on scales A divided by the rest en- 
ergy density of the shock upstream. Note that e# depends on 
the length scale A and is not directly related to the confining 
magnetic field energy density measured on scales Ad e f • 

We define x = A/(R/ 8T 2 ) so that the unperturbed motion of 
an upstream fluid element toward the shock follows the char- 
acteristic dx/dr = -c/(R/8T 2 ). Ignoring the amplification of 
the magnetic field during the linear phase of the nonthermal 
proton-driven streaming instability (§ I2.41 i we obtain 



d 2 i 
dx 2 



0.02 fiee^E^egt 
r4/3^i/6 m ^ c 7/3 >maxA - 



(15) 



Given that the upstream fluid starts from rest, £(1) = 
d^/dx{\) = 0, equation (TT3T > can be integrated to find the trans- 
verse displacement 



0.02 /Lte4 /2 gto{ 3 e n t 

r4/3^1/6 m f c 7/3 7pmaxA 



g(x)sin0. (16) 



where g(x) = xlnx-x + 1 for the coherent field and g(x) = 
\x^l 2 -2x+ § for the tangled field. Setting fj, = 1 and A = 10 
we obtain 



R 



sr v 1.6 x 10 4 



-1 . , , 1/2 / p \ 1/3 



(— V 

Vio-v 



100 



1/3 



1 cm 



V10 53 erg 

1/6 

, ; gCxOsin^ 



(17) 



Thus for 7 Piinax ^> 10 4 , the distance traveled by the upstream 
fluid is <C R/8T, which is the width of the blastwave in the 
frame of the shock downstream. The value of £(x = 0) de- 
creases with increasing 7 z ,. max because when protons are ac- 
celerated to higher energies, the number of nonthermal pro- 
tons between which the available energy ^ e n t£tot is parti- 
tioned decreases, while the lower energy nonthermal protons 
are confined closer to the shock, both of which decrease the 
density of nonthermal protons far from the shock where the 
Ampere force-driven acceleration takes place. 

Under what circumstances will the transverse displacement 
result in a nonlinear density contrast on scales A? This re- 
quires £(x) ~ A before the fluid element is overtaken by the 
shock at x = 0. However, before the fluid arrives at the shock 
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transition, some assumptions entering the derivation of equa- 
tion ( fTo*b may be compromised: its motion may become rel- 
ativistic (d^/dr ~ c), it may enter the region in which the 
pressure transfered by the nonthermal protons to the upstream 
fluid starts to accelerate the upstream fluid in the original rest 
frame of the shock upstream, and it may also enter the region 
containing accelerated nonthermal electrons where the return 
current vanishes. The latter problem arises when x ~ x coo i 
where x coo i = A c (7^ max )/(,R/8r 2 ) is given in equation (f2j). It 
is evident from equation OJ that the return current vanishes 
only in the region populated by the lowest-energy nonthermal 
protons (7 P < 10 3 ), which for 7 p , max ^> 10 3 , as expected, will 
be confined very close to the shock transition. We are not able 
determine the value of x at which the upstream starts to be ac- 
celerated in its own frame but we do assume that it is smaller 
thanx coo i. 

To assess the distance at which the transverse motion of 
the fluid becomes relativistic, we focus on the special case 
of a tangled deflecting field; the case of a coherent deflect- 
ing field is qualitatively similar. The time at which the 
upstream motion becomes relativistic is obtained by setting 
dijdx = -R/8T 2 (i.e., d£/dr = c) to find 

Xml= (l-lEj^] (tangled field) , (18) 



7rel,crit , 

for 

7p,max < 7rei,criti where 

0.4 T^W^E^t 



7rel.crit ' 



cV3 m y 6 A 



-3 x 10 b 



100 



10-V \ 10 53 erg/ V0.1 



1/3 



\ 4/3 / n y/6 
J I Tot? J 



(tangled field) . (19) 



When 7p. m ax > 7rei,crit> me fluid does not attain relativistic 
transverse velocity before it reaches the shock. Then, the 
numerical estimate in equation ( fTTI i at x = where g(x) ~ 1, 
which was derived assuming Newtonian motion, shows that 
the transverse displacement decreases with 7 ; ,. max . 

Conversely, for 7 /) , max <C 7rei,crit> the fluid motion formally 
becomes relativistic well before the fluid reaches the shock. 
The maximum transverse displacement that can be attained in 
the Newtonian regime can be evaluated by substituting x = x K \ 
in equation (fT~6b to obtain 



£(x re i)~0.3 



7w 



jiee^g 2 n} l 2 T A sin^> 

0.004 -f ^) (-2L) 
8T V 10 6 / V10-V 

emV 1 / gtot 
0.1/ V!0 53 erg 



-1/2 



r 
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-7/3 



1 cm- 



-1/6 



(tangled field) 



(20) 



for 7 P ,max *C 7rei,crit- Acceleration of the fluid may not stop in 
the Newtonian regime; indeed, it will likely continue into the 
transrelativistic regime, where the equation of motion is sim- 
ply d^/dr ~ c. Then, however, the motion will not be purely 
perpendicular to the direction of shock propagation and the 
fluid will start accelerate, in the upstream frame, in the par- 
allel direction. Here we do not take into account this behav- 
ior and conservatively restrict our attention to the maximum 



transverse displacement that the fluid reaches while its motion 
remains Newtonian. 

Since the transverse displacement increases with 7 Pimax 
when the transverse motion becomes transrelativistic and de- 
creases with 7 /J , max when it does not, the largest Newtonian 
displacement is possible when the maximum proton Lorentz 
factor is about equal to 7 re i,crit in equation ( fT9l and is simply 
given by £ max ~ ^R/8T 2 , which is much larger than the pro- 
ton plasma skin depth c/uj p = c/{A%e 2 n/m^) 1 ' 2 of the shock 
upstream. Of course, the maximum transverse displacement 
is achieved only if a magnetic field that reverses on scales 
A ~ Cmax is present. We speculate about the origin of the re- 
versing field in § 12.41 below. Note that £ max is close to the 
maximal value allowed by causality, R/8T 2 . 

2.3. Density Inhomogeneity, Vorticity, and Dynamo 

Here, we study the evolution of an initially uniform density 
fluid in response to the Lorentz force-driven transverse mo- 
tion discussed in §12.21 If the magnetic field reverses itself on 
scales ~ A, the average linear density contrast arising from 
the transverse motion, for linear displacements, («A, will 
be given by 



P-P 



(21) 



We expect nonlinear density contrast on scales A provided 
that £(x max ) > A, where x max = max{x coo i,x re i}. This condi- 
tion places an upper limit on the magnetic field reversal length 
scales A on which a nonlinear density contrast can build up 
before the fluid reaches the shock. 

When fluid with nonlinear density contrast reaches the 
shock transition, provided that the width of the density jump 
of the shock transition is smaller than the length scale asso- 
ciate d with the contrast, vortic ity is generated at the transi- 
tion. ISironi & Goodmanl d2007l) have calculated the fraction 
of shock energy that is converted into vortical energy for an 
ultrarelativistic blastwave propagating into a clumpy medium. 
They provide a fitting formula for vortical energy fraction e vort 
as a function of clump size L, peak density contrast (5 max , and 
volume-filling factor n c iump£ 3 , which reads 

'0.36« c i ump L 3 \ 5l 



5 2 

max 



1+0.176 (<W)i 



(22) 



For L ~ A and n c i ump ~ A 3 and mildly nonlinear density con- 
trast <5 max ~ 1 , this implies e vort ~ (3T)" 1 . Th e collisions of and 
compression in accelerated plasma shells (Bell 2004, [2005b 
may amplify the density contrast into the strongly nonlinear 
regime <5 max ^> 1, but we do not attempt to estimate the mag- 
nitude of this effect. 

If vorticity is generated on scales A <C R/T, the result- 
ing eddies can amplify the downstream magnetic field to 
levels Cb down ~ e vnrt v i a a turbulent dynamo mechanism 
dMeneguzzi et all 119811; ISchekochihin et al.l 12004 . Strong 
amplification is expected when the number of vortical eddy 
turnovers in the shock downstream is larger than unity, but this 
is automatically expected when the density contr ast i s nonlin- 
ear on scales <C R/8T. Indeed, our analysis in § !2.2l suggests 
that nonlinear density contrast will be generated on scales 
A < £ max ~ |7?/8r 2 , for which the expect ed number of vorti- 



cal ed dy turnovers is 3> 3T (see § 4.2 in ISironi & Goodmanl 
120071) . On the other hand, since £ max is much larger than 
the plasma skin depth c/u> v , the magnetic field that the tur- 
bulent dynamo generates will not be susceptible to the fast 
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collisionless decay discussed by Chan g^tal] (120081) that af- 
fects the much smaller scale field that kinetic instabilities can 
generate in the shock t r ansition layer (e .g. 



Frederiksen et alj 12004 iMedvedev et all 120051: iKeshet et al. 
2008L ISpitkovskyll2008bl and references therein). 



2.4. Reversing Field from Bell's Streaming Instability 

In § 12.21 we assu med that the mag netic field reverses di- 
rection on scales A. Bell ( 20041 120051) has shown that an in- 
stability driven by the nonthermal ion streaming can convert 
an initially unif orm field into a reversing, helical magnetic 



field (see, also. iReville et all l2006HZirakashvili et al.l 12008b 
Pelletier, Lemoine, & Marcowith 2008). For ultrarelativisic 
streaming and helical perturbations of the form £, + i£ y oc 
exp i(kR - cot) (x and y are two coordinates perpendicular to 
R), the instability has dispersion relation 



2 t kJretBR 

-k v A 

pc 



(-£)- 



o. 



(23) 



where va is the Alfven velocity and Br = B • R. The fastest 
growing mode has wavelength and growth rate 



A 



Br 



fast : 



Tf as 



Mitp 



fast — en ntp 



/ nm p 

\ 7T 



-1/2 



(24) 



respectively. 4 For nonlinear density contrast production on 
scales Af as t we need £(x coo i) > Af ast and £(je re i) > Af ast as well 
as Tf ast > %Y 2 cjR. However, the wavelength of the fastest 
growing mode is very small if B is the microgauss field typical 
of the interstellar medium. All longer wavelengths A > A tast 
also grow on rates 



T(A) - T 



fast 



A 

Afast 



-1/2 



(25) 



This seems to suggest that a reversing field may be generated 
on a wide range of scales. 

The instability has enough time to reach nonlinear growth 
only when the growth rate exceeds the inverse crossing time 
of the shock precursor 



T(A)» 



R 

8T\ 



(26) 



To evaluate the growth rate, we substitute equations ( fTOl co- 
herent field) with A = R/ST 2 and d24ji in equation (|25]l to ob- 
tain 

2.4 e V2 c VV/4 e y4 e V4 r 2 

T(A) J,: . (27) 

7;^axAV2 m y 4 A'/2 

Where eg parametrizes the strength of the pre-existing coher- 
ent field, not to be confused with the same symbol expressing 
strength of the length scale-dependent tangled field in previ- 
ous sections. The condition for nonlinear growth in equation 
is then fulfilled when 



7w 



8r 2 \4.6x 10 6 

1/6 



1M( J_V /3 



(— 
V 1 cm 



(. 



10- 9 / vo.iv \100J 

1/3 



cm 3 / \ 10 53 erg 



(28) 



4 It is a peculiar characteristic of ultrarelativistic blastwaves that the gy- 
roresonant wavelength is always longer than the particle's flight path in the 
shock upstream. 



which suggests that nonlinear growth is very broadly expected 

On Scales A ^ £max ~ R/ET 2 When 7/?.max ^5 7rel,ciit- 

When 

7p,max ^S> 7rei,crit> nonlinear growth is still expected on scales 
smaller than R/8T 2 but still much larger than the plasma skin 
depth. 

3. DISCUSSION 

The complexity of the structure of ultrarelativistic collision- 
less shocks stems from the multiscale nature of plasma self- 
organization. This and the preceding analytical and numer- 
ical analysis have identified organization on various spatial 
scales between the proton plasma skin c/lu p ~2x 10 7 cm and 
the width of the blastwave in the upstream frame R/8T 2 ~ 
10 12 - 10 18 cm. While the dynamical organization of the 
plasma on different scales is interdependent, our understand- 
ing of the scale interdependence in the precursors of ultrarel- 
ativistic collisionless shocks is still in its infancy. Motivated 
by the present and other preliminary investigations (see, e.g., 
Milosavlie yic & Nakari 120061: iKatz et all 120071: IKeshet et ail 
120081: ISpitkovsk v 2008b), we propose several possible regions 
within the structure of an ultrarelativistic shock precursor in 
which the physical processes that drive plasma organization 
change character as a function of distance from the shock. 
We, however, caution that not all of these regions have to be 
present in all — or perhaps any — of the blastwave precursors. 
Similarly, we are not able to determine the relative ordering 
of the regions, with respect to their distance from the shock 
transition, with any degree of certainty. In a shock that pro- 
duces a power-law spectrum of accelerated nonthermal parti- 
cles (see iUfor potential limitations of this key assumption), 
we propose that the shock precursor may contain the follow- 
ing concentric regions: 

1 . An outermost region extending to the maximum possi- 
ble radial distance from the shock, measured in the reference 
frame of the shock upstream, of ~ /?/8T 2 , that is causally as- 
sociated with the explosion. This region contains a low den- 
sity of nonthermal protons in the high-energy tail of the non- 
thermal proton spectrum. Because inverse Compton cooling 
losses prohibit the existence of nonthermal electrons in this 
region, a return current neutralizing the charge separation cre- 
ated by the high-energy accelerated nonthermal protons flows 
in the upstream plasma (§ 12.11 ). The effect of the return current 
is twofold: (a) starting from an initial, weak, coherent, quasi- 
uniform magnetic field, the return current nonresonantly ex- 
cites circularly-polarized Alfven waves, thereby converting 
the field component parallel to the direction of shock propaga- 
tion into a wound, perpendicular, helical component (§ 12.41 ); 
(b) given a perpendicular magnetic field, which could be pre- 
existing or generated by the return current, the Ampere force 
associated with the return current acce lerat es the upstream 
plasma in the perpendicular dire ction (§ 12.21 ) and produces a 
density contrast in the plasma (§ 12.31 ). More speculatively, it 
has been hypothesized that besides generating nonlinear den- 
sity contrast, the return current interaction significantly am- 
plifies the pre-existing magnetic field; this possibility awaits 
verification in numerical simulations. Because the density of 
nonthermal particles in this region is low, the region is either 
stable to small-scale plasma instabilities, such asjthe trans- 
verse Weibel (or filamentation) instability (e.g., lFriedl[T959h 
and its relatives, or more likely, the small-scale instabilities 
are present but their growth times are longer than the crossing 
time of the region < R/8T 2 c. 5 This region should exist unless 

5 Small scale instabilities may also occur of the upstream medium is en- 
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the compression layer, which is the region 4 below, extends 
all the way to the edge of the maximum radial distance from 
the shock (~ R/8T 2 ), which renders the shock transition layer 
as wide as the downstream. 

2. A region still devoid of nonthermal electrons and en- 
dowed with the return current, but where the density of non- 
thermal protons is high enough to drive the small-scale plasma 
instabilities, with wave vectors k^oj p /c. At present, the dom- 
inant unstable kinetic mode in the regime in which a very low 
density beam of particles with an ultrarelativistic forward mo- 
mentum and an ultrarelativistic transverse momentum spread 
streams through a cold background plasma has not been iden- 
tified. It is possible that the small-scale instabilities can gen- 
erate an electromagnetic field efficiently and that the field 
enhances the perpendicular velocity diffusion of nonthermal 
particles, which is a necessary element of the DSA cycle. The 
magnetic field in this region is a combination of a small-scale 
field generated by the instabilities and a larger scale field gen- 
erated from the pre-existing field by the return current. If 
the latter is strong enough, depending on its resulting orien- 
tation relative to the shock normal, it may suppress some of 
the small-scale instabilities, or it may drive small-scale turbu- 
lence in the flow, especial ly if the resulting field is per pendic- 
ular to the shock normal dHededal & Nis hikawa 2005). 

3. A region close enough to the shock transition that the 
inverse Compton cooling time of the nonthermal electrons is 
longer than the duration of the DSA cycle. Since this region is 
populated by nonthermal electrons and protons, the nonther- 
mal particle precursor is neutral overall, and no return current 
is present. However, small scale instabilities are still present 
and are more intense than in the preceding regions. The small 
scale instabilities generate a strong magnetic field in the up- 
stream plasma, which, when advected to the shock transi- 
tion, facilitates upstream particle scattering and thermaliza- 
tion. The small scale field may also play a critical role in the 
injecti on and acceleration of particles in the DSA mechanism 
(e.g jLemoine et alX 2006). The relative ordering of regions 2 
and 3 may vary. Depending on the parameters and initial con- 
ditions, it may happen that the region containing nonthermal 
electrons extends farther from the shock than that in which 
small scale plasma instabilities excited by nonthermal parti- 
cles play an important role in setting the shock structure. 

4. The shock transition region, where the upstream plasma 
is compressed. At the leading edge of this transition layer, 
the compression may result from an effi cient momentum 
transfer from the no nthermal particles (e.g., Bla ndfordll 19801 : 
iDrury &"Vo lk 1981), which could be coupled to the upstream 
plasma by the small scale instabilities. In this region, the up- 
stream plasma starts to accelerate (in the rest frame of the 
far upstream) and decelerate (in the rest frame of the shock). 
The final stage of compression, where the upstream becomes 
Newtonian in the rest frame of the shock, takes place when the 
upstream interpenetrates the downstream isotropized plasma, 
and so the hydrodynamic jump conditions are realized. If the 
nonthermal particles play a role in initiating the compression, 
than the width of this layer is much larger than the plasma 
skin depth. If, however, efficient compression commences 
only where the upstream interpenetrates the downstream ther- 
malized plasma, the width of the compression layer may be 
as small as of the order of the skin depth. While in the 
preceding regions we expect a decrease in the mean mag- 
netic field correlation length (A) = (Afi 2 (A)) / (B 2 (X)} as the 

riched by an -pair cascade jRamirez- Ruiz et alJ2007l) . 



plasma approaches the shock, in the shock transition region 
we expect an opposite trend: the small scale instabilities have 
reached saturation and the small field they have generated de- 
cays via a kinetic da mping mechanism as the shocked plasma 
travels downstrea m (Gruzinov 2001; Medvedev et al.l 12005b 
IChang et al.ll2008UKeshet et al]l2008brSpitkovskyl2008bl) . As 
the small scale field decays, the residual field, which par- 
takes jrrthejiyjichroto^ (see, 
e.g., iPanaitescu & Kumarl 12001 1 Yost et al.l 120031) and must 
be much stronger than the shock compressed pre-existing 
mangetic field of the normal interstellar medium (see, e.g., 
Gruzi novll2001l: lPiranll200 5b). is the larger scale field gen- 
erated in regions 1-3 above. The shock transition region is 
also the site of the onset of the e stablishment of electron-ion 
equipartition dSpitkovskvl2008a]) . 

Finally we note that if the highest energy nonthermal parti- 
cles traverse regions 1-3 during a DSA cycle, that is, these 
particles are deflected by the magnetic field of the non- 
acc elerated upstream, then an upper limit on j p . lwdx can be 
set (Milosavljevic & Nakar 2006). This limit is derived from 
the constraint of causality on the self organization scale of 
the upstream magnetic field (~ R/8T 2 ) and the maximal 
magnetic field energy density in the non accelerated up- 
stream (~ nm p c 2 ). This upper limit is, roughly, 7 p , max < 
10 10 n 1/2 (/?/10 18 cm) implying the GRB external shocks 
are not promising sources of ultra-high energy cosmic rays 
(UHECR) with Lorentz factors, under the assumption that 
they are protons, of ^ p > 10 11 . 

4. CONCLUSIONS 

We have examined the influence of a pre-existing magnetic 
field on the structure of ultrarelativistic shock waves, such as 
in the external GRB shocks that may be behind the GRB after- 
glow phenomenon. From the hypothesis that the shock wave 
accelerates charged particles to high Lorentz factors, we were 
able to derive the following conclusions: 

In GRB external shocks that accelerate ions to Lorentz fac- 
tors 7, ^> 10 3 , inverse Compton cooling in the shock up- 
stream hinders nonthermal electron acceleration to high ener- 
gies. Because of this, nonthermal electrons do not penetrate as 
far into the shock upstream as the highest-energy nonthermal 
ions do. Thus, a region exists in the shock upstream far from 
the shock transition at R ~ ct[l -(8r 2 ) -1 ], but still behind the 
photon shell at R = ct, where return current flows to cancel the 
charge separation between the nonthermal electrons and ions. 

As the return current flows across the weak pre-existing 
magnetic field of the shock upstream, the Ampere force ac- 
celerates the upstream plasma perpendicular to the direction 
of shock propagation. If the upstream field reverses on scales 
A <§C R/8T 2 , the transverse acceleration produces a nonlinear 
density contrast on the same scales. If the pre-existing field 
is quasi-uniform (non-reversing), Bell's streaming instability 
will efficiently convert the quasi-uniform field into a reversing 
field during the light-crossing time of the shock precursor. 

As the shock transition sweeps the resulting density clumps, 
the bulk kinetic energy of the blastwave is converted into vor- 
tical energy. The vortical eddies in the shock downstream 
turn over many times and amplify the magnetic field via a 
dynamo mechanism. For clumps of peak density contras t 
of about one, the vortical energy ( Sironi & Goodman 2007) 
and with it the magnetic energy density can be amplified to 
e B ~ (3TT\ where in GRB external shocks, T < few x 100. 
This is consistent with the observed GRB light curves and 
spectra. The generated magnetic field will be coherent on 
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scales much larger than the plasma skin depth and will, as 
such, not be susceptible to rapid collisionless decay. Taken 
together, these results indicate that the microgauss magnetic 
field (eg ~ 10~ 9 ) typical of the interstellar medium is suffi- 
cient to drive Bell's mechanism in ultrarelativistic GRB ex- 
ternal shocks and deliver nonlinear density inhomogeneities 
to the shock transition. Therefore, magnetic field amplifica- 
tion through shock vorticity generation is expected even when 
the circumburst medium is not initially clumpy, and we expect 
it to operate in long and short GRBs. 

We speculatively propose a structure for the shock pre- 
cursor consisting of four concentric zones characterized by 
different forms of plasma self-organization. First, the out- 
ermost region of the precursor is the region populated by 
streaming, nonthermal protons but devoid of nonthermal elec- 
trons. There, the plasma develops non-linear density inhomo- 
geneities. Second, a region that still contains a return cur- 
rent but where the nonthermal particle precursor streaming 



through the shock upstream is susceptible to small-scale ki- 
netic plasma instabilities. Third, a region containing the rela- 
tively low-energy nonthermal electrons that can be acceler- 
ated in the shock in spite of the inverse Compton cooling 
losses; in this region, the return current is absent. Finally, 
the fourth region is the shock transition itself, where kinetic 
instabilities and plasma waves facilitate plasma isotropization 
and establish the hydrodynamic jump. The relative ordering 
of the second and the third region may vary. 
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